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Abstract. Hubble Space Telescope (HST) limits on the 
amount of halo dark matter (DM) in the form of very 
low-mass (VLM) stars close to the hydrogen-burning limit 
have been interpreted as excluding such stars as viable 
candidates. However, these limits assume that the stars 
are smoothly distributed and have at least the metallicity 
of Population II stars, whilst some baryonic DM forma- 
tion theories predict that they may instead be clumped 
into globular-cluster configurations and have close to zero 
metallicity. I re-analyse the HST data employing the zero- 
metallicity VLM star models of Saumon et al. (1994), 
which predict V — I colours below the cuts of previous 
analyses for stars below 0.2 Mq. From the models I de- 
rive new limits on the allowed halo fraction comprising 
VLM stars for both the unclustered and clustered cases. 
In the unclustered regime I find a 95% confidence up- 
per limit on the allowed halo fraction of 1.4% inferred 
from 20 HST fields, comparable to limits derived by pre- 
vious studies for non-zero metallicity populations. In the 
cluster scenario I show that clusters of mass M and ra- 
dius R can satisfy both HST and the recent MACHO 
gravitational microlensing results, which indicate a lens 
halo fraction of 40% for a standard halo model, provided 
R 5s 1.2(M/10 4 Mq) 74 pc. However, existing dynami- 
cal limits restrict the allowed range to a tiny region char- 
acterised by M ~ 4 X 10 4 M and R ~ 3 pc. Further- 
more, consistency between MACHO and HST demands a 
present-day clustering efficiency of 92% or better. Intrigu- 
ingly however, the cluster mass implied by these limits is 
theoretically well motivated and the VLM star scenario 
may also help to provide an explanation for the faint red 
'halo' light recently reported around another galaxy. 

Key words: stars: low-mass, brown dwarfs - globular 
clusters: general - Galaxy: halo - dark matter 



1. Introduction 

Recent advances in attempts to detect or constrain the 
nature of dark matter (DM) in the halo of our Galaxy have 
led to a variety of constraints on both baryonic and non- 
baryonic candidates (Carr 1994; Jungman et al. 1996). 

The microlensing r ate o bserved by the MACHO collab- 
oration (Alcock et al. 1997 ) towards the Large Magellanic 
Cloud (LMC) indicates that the likely halo DM fraction 
in compact form is around 40%, and comprises objects 
in the mass range 0.1 — 1 Mq, though the inferred frac- 
tion and mass range are sensitive to the assumed halo dis- 
tribution function. This therefore suggests that the halo 
may comprise a roughly equal mixture of baryonic and 
non-baryonic matter, though there are non-baryonic can- 
didates, such as primordial black holes and shadow mat- 
ter, which can also explain the microlensing events. In any 
case, the MACHO results imply that no single candidate, 
baryonic or non-baryonic, can explain all of the halo DM 
unless the halo distribution function departs significantly 
from the usual assumption of an isothermal sphere. 

If the DM responsible for the observed microlensing 
is baryonic then the inferred mass range implicates cither 
white-dwarf remnants of an early generation of stars or 
very low-mass (VLM) stars close to the hydrogen-burning 
limit. However, in our Galaxy the number density of halo 
white dwarfs is strongly constrained by the present-day 
helium and metal abundances of the interstellar medium 
(Carr et al. 



1984; Ryu et al. 1990; Adams & Laughlin 



1996) and, assuming the white dwarfs are younger than 
18 Gyr, by number counts of high- velocity white dwarfs 
(Chabrier et al. 199£). Counts of high-redshift galaxies 



also appear to rule out white dwarfs from contributing 
significantly to halo DM in other galaxies (Chariot & Silk 



1995) 
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VLM stars are currently even less favoured since pho- 
tometric surveys appear to place even st ronge r limits on 
their contri bution (e.g. Gilmore & Hewett 1983 Richstone 
et al. |1992[ Bahcall et al. |l994| , hereafte r BFG K; Graff & 
Freese |99§ hereafter GF; Flynn et al. |l996[ ). The latest 
MACHO results also appear to argue against most of the 
DM comprising sub hydrogen-burning brown dwarfs un- 
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less the DM is distributed more like a maximal disc, in 
which case any halo must have a very low total mass inte- 
rior to the LMC in order to remain consistent with local 
surface-density and rotation-curve constraints. Within the 
context of standard halo models therefore, there is great 
difficulty in explaining the high halo fraction detected by 
MACHO by invoking a single baryonic candidate. 

In this paper I re-examine the VLM star hypothesis 
for the case where the stars have zero-metallicity, as may 
be the case if the halo is formed from primordial matter. 
In this case their colours will be different than assumed 
by previous studies, so constraints from these studies need 
not necessarily apply to zero-metallicity stars. I also exam- 
ine the possibility of DM clustering, and how that effects 
number counts and microlensing statistics. Such clustering 
is predicted by some DM formation theories. 

The plan of the paper is as follows: Sect. || calculates 
the space density constraints for a smoothly distributed 
population of zero-metallicity VLM stars imposed by num- 
ber counts from 20 of the 22 HST fields obtained by 
BFGK. These calculations adopt the latest numerical pre- 
dictions for zero-metallicity VLM star photometry. It will 
be shown that the limits on such a population are compa- 
rable to those for populations of non-zero metallicity, with 
the allowed halo fraction being at most 1.4%. Section || ex- 
tends the analysis to include the effect of clustering and, 
using a combination of the HST number-count limits, MA- 
CHO observations and dynamical constraints on the per- 
mitted mass and radius of clusters, places upper limits on 
the cluster space density. It will be shown that whilst dy- 
namical constraints strongly limit the allowed mass and 
radius of clusters, they nonetheless permit a cluster DM 
fraction which can simultaneously explain the MACHO 
and HST observations, though both the required mass M 
and radius R need to be finely tuned (M ~ 4 x 10 4 Mq, 
R ~ 3 pc) . I further show that the microlensing signatures 
of these clusters are essentially indistinguishable from the 
unclustered scenario, and in particular are not expected to 
modify the inferred halo fraction. Finally in this section, I 
use the halo DM fraction inferred by MACHO to compute 
the efficiency with which VLM stars must be clustered in 
order to remain compatible with HST observations. The 
required efficiency, better than 92% at the present day, 
places the strongest demand on the scenario. The main 
findings of the paper are discussed in Sect. ^. 

2. Limits on zero-metallicity VLM stars 



Saumon et al. (1994) have computed a grid of fully non- 
gray atmospheric models for zero-metallicity (Z — 0) 
VLM stars and brown dwarfs in the mass range 0.01 — 
0.2 Mq. Their calculations take account of all the ma- 
jor sources of opacity for a pure H+He mixture includ- 
ing H2, H and H~, a s well as other sources tabulated by 
Lenzuni et al. ( 1991 ). The resulting spectral energy dis- 
tributions are found to deviate significantly from black- 



body below an effective temperature of 4000 K. Saumon 
et al. also present predictions for the photometry of VLM 
stars (0.092-0.2 Mq) from M through to V bands which 
show that these stars would appear significantly bluer 
than stars with metallicity comparable to that measured 
for the Galactic spheroid (Z ~ 0.01 Zq) or disc (Z ~ Zq) 
populations. 

In particular, Saumon et al. find that the expected 
V— I colour for zero-metallicity VLM stars ranges between 
1.27 for 0.2-Mq stars to 1.57 at the hydrogen-burning 
limit (0.092 Mq). This is somewhat below the threshold 
V — I values assumed by BFGK, GF and Flynn et al. 



( 1996 ) in their analyses, whose inferences arc based on 
the number counts of stars with V — I ib 2. Whilst one 
can not say for sure whether any dark halo population has 
a metallicity substantially less than that of the spheroid, 
there is clearly a need to re-examine the HST data using 
the zero-metallicity predictions of Saumon et al. in order 
to close up this last remaining 'loophole'. 

The positi ons o f the 22 HST fields are tabulated in 
Gould et al. ( |1996[ ), along with the corresponding min- 
imum and maximum /-band magnitudes for each field. 
(All selections were performed in the / band.) Note that 
two locations (I = 82°, b = -19° and I = 134°, b = -65°) 
each have 2 fields in very close proximity. 

The limiting maximum /-band magnitude, which de- 
termines the maximum distance out to which a VLM star 
will be seen, ranges from 22.56 to 24.40 and is deter- 
mined for each field according to the ability to discrim- 
inate clearly between stellar and extended sources. The 
minimum /-band magnitude ranges from 17.05 to 19.45 
and sets the minimum distance at which the stars can be 
satisfactorily imaged. Objects with magnitudes below this 
limit produce saturated images. It should be noted that 
BFGK calculate the /-band magnitude limits from the 
HST /' band (F814W filter) assuming the stellar spec- 
tral energy distributions listed by Gunn & Stryker ( 1983 ). 
They therefore do not strictly apply to Z — VLM stars, 
though any differences in calibration will be small and are 
therefore neglected here. 

In computing the minimum and maximum observable 
distances for each VLM star mass, and for each field, I 



follow Gould et al. (1996) in converting the B-band ex- 
tinction values determined by Burstein & Heiles ( [1982 ) to 
/-band reddenings. The extinction is assumed to be con- 
fined to a disc of thickness \z\ = 100 pc. Within this disc 
the extinction per unit distance is taken to be uniform. 

The halo mass within the volume defined by the min- 
imum and maximum distances d m in, tax, and the solid 
angle per field O = 3.7 x 10~ 7 sr, is calculated assuming 
a spherically-symmetric softened isothermal halo density 
distribution of the form 



p(x,l,b) = po 



Rl 



2xRq cos I cos b + a 2 



(i) 
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where x is the distance measured along the observer's line 
of sight, / and b are Galactic coordinates, Rq = 8 kpc 
is the Sun's Galactocentric distance, po = 0.01 Mq pc~ 3 
is taken to be the local DM density normalisation and 
a = 5 kpc is the assumed halo core radius. For the small 
solid angles considered here, this gives an integrated halo 
mass between d m -, n and d mFn < for field i of 
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x 2 p(x, li, bi) dx. 



(2) 



The integral can be performed analytically though the re- 
sulting expression is long. 

Note that the integral limits in Eq. |2| are implicit func- 
tions of the VLM star mass m. If one assumes that VLM 
stars have masses to* > to then a lower limit on the ex- 
pected number of detectable VLM stars is 



VLM 
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(> TO) > — y*M ha lo,iO 
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d m in,i ( m ) 



x 2 p(x, bi) dx 



(3) 



for n independent fields. The direction of the inequality 
reflects the fact that the dependency of M hal0ji on to is 
steeper than the first power of to. Fields 1 and 19, using 
the order in which they are listed in Tab. 1 of Gould et al. 

are discarded in this analysis because of their close 
proximity to fields 2 and 20, respectively. Whilst fields 1 
and 2 do not actually overlap, field 1 is nonetheless ex- 
cluded here to provide consistency with the cluster anal- 
ysis of Sect. ^2, where statistical independency of neigh- 
bouring fields is an important criteria. ./Vvlm is therefore 
summed over n — 20 rather than 22 fields. 

The median value for c? m j n is found to range from 
330 pc for O.O92-M stars up to 1 kpc for 0.2-Mq stars. 
For ci max the median values are 3.3 kpc and 10.1 kpc for 
0.092-Mq and 0.2-Mq stars, respectively. 

Applying Eq. 3 to VLM stars at the hydrogen-burning 
limit mass of 0.092 Mq, one finds an expectation num- 
ber of stars in the HST fields of 6 310, and for O.2-M 
objects the expectation is nearly an order of magnitude 
larger at 60 100 stars. These numbers take account of the 
fact that data from one-third of field 4 (i.e. data from one 
of the three detector chips) had to be discarded by Gould 
et al. ( 1996j ) due to problems with receiving the data from 
HST, and that as much as 2% of each of the fields was 
discarded due to emission from background galaxies. The 
number of stars detected in the 20 fields with V — I val- 
ues in the range 1. 2 to 1.7, spanning the range predicted 
by Saumon et al. ( 1994 ), is only 75. The 95% confidence 
level (CL) upper limit on the average, for a realisation of 
75 stars, is 91. Therefore, even if one assumes that all of 
the stars detected by HST are halo VLM stars right on 
the hydrogen-burning limit, their contribution to the halo 
DM can be no more than 1.4% at the 95% CL, and the 



limit is correspondingly stronger than this for more mas- 
sive objects. In fact it is likely that a significant fraction 
of these objects may belong to the disc or spheroid. 

The limit of 1.4% is stronger than that inferred by 
BFGK for solar-metallicity VLM stars from their anal- 
ysis of one of the HST fields. Their results translate to 
a 95% CL upper limit of less than 4% for stars with 
V - I > 3, adopting p = 0.01 M pc~ 3 . The 95% up- 
per limit inferred by GF for low-metallicity stars with 
2 < V — I < 3 corresponds to a halo fraction of less 
than 0.9% for the same field. However, both of these stud- 
ies exclude from their analyses stars with V — I < 2, so 
the limits presented here, which are derived from 20 fields, 
constitute a completely independent check on previous re- 
sults. The conclusion is that, regardless of their metallic- 
ity, VLM stars do not contribute significantly to the halo 
DM, at least under the assumption that they are smoothly 
distributed in the halo. 



3. VLM star clusters 

3.1. Rationale 

The rationale for invoking clustered DM comes from a 
number of considerations. Firstly, visible globular clus- 
ters are observed out to great distances from the Galac- 
tic centre suggesting they may have an important role 
in halo formation, though they are identified with the 
spheroid population. Secondly, some of the most promis- 
ing theories for the formation of l arge amounts of com- 



pact halo baryonic DM (Ashman 199C ; De Paolis et al 



1995 ) predict that the DM should be clustered into dark 
globular c luste rs with a typical mass of around 10 4 Mq 
(Ashman 1990 ). These theories predict that the clusters 
should comprise either brown dwarfs or VLM stars, due 
to high gas pressures suppressing the minimum fragmen- 
tation mass. They also explain the typical mass and spa- 
tial distribution of the visible cluster population since 
the theories are extensions of the Fall-Recs theory of 
globular-cluster formation (Fall & Rees 1985). The baryon 
Jeans mass at the cosmological epoch of recombination is 



Mi 



BJ 



1.3 x 



lo^B^yv 



Mq, which is also close to 



10 4 Mq for a cosmological density Slo of order unity and a 
baryon density fie satisfying cosmological nucleosynthetic 
constraints: 0.01 < fish 2 < 0.024, where h is the Hub- 



ble constant in units of 100 km s^ 1 Mpc^ 1 (Krauss 1995). 
Lastly, limits on the space density of clusters provide firm 
upper limits on the allowed space density of VLM stars, 
provided one assumes that the VLM star distribution fol- 
lows the halo density on the large-scale average (i.e. scales 
larger than the typical cluster separation). 

The effect that clustering can have on number-count 
statistics is well illustrated by considering the expectation 
distances both to the nearest unclustered VLM star and 
to the nearest cluster. Let us assume that VLM stars con- 
tribute some fraction /h to the halo DM and that their 
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density is everywhere fhp, where p is given by Eq. §. Let 
us also take the fraction of VLM stars residing in clusters 
to be everywhere f c , so that the fraction of unclustered 
VLM stars is correspondingly 1 — / c . Then, for VLM stars 
with mass m, the local unclustered number density is 



/h(l - /c)po 



■m 



o.i/ h (i-/ c 



0.1 M 



pc 



(4) 



whilst the local number density of VLM star clusters of 
mass M is 



n c ,o 



/h/cPO 

M 

= io- 6 A/c 



M 



10 4 M 



© 



pc 



(5) 



The expectation value for the distance d to the nearest ob- 
ject appearing within a survey of solid angle O, a ssum ing 
a local number density no, is (c.f. Kerins & Carr 1994) 



(d) = Qrio / x 3 exp(— QriQX 3 /3) dx 
Jo 

= (3/0) 1 / 3 L(4/3)< 1/3 , 



(6) 



where L(4/3) = 0.893 is the gamma function. So, using 
Eq. 4 and Eq. 5, the expectation distances to the nearest 
unclustered VLM star and star cluster are, respectively, 



(du) = 



0.14/- 1/3 (l-/ c )- 1 /3 

n 



{d c ) = 6.6/^7, 



7.4 x lO- 6 sr 

l/3 f -l/3 



1/3 / x 1/3 

m x 



0.1 M r . 







kpc (7) 



O 



7.4 x 10" 6 sr 



-1/3 



M 



10 4 Ma 







1/3 



kpc, (8) 



where Q is normalised to a solid angle 20 times that of a 
single HST field. 

The large value for (d c ) indicates the need to consider 
the variation of p along the line of sight in order to obtain 
a more accurate answer, though the estimate of Eq. 8 is 
sufficient to show that one might not see any clusters at 
all within the HST fields, even if they comprise all the 
halo DM. Additionally, the expected fluctuation in source 
counts from field to field will be governed by the number 
density of clusters and not the spatially-averaged density 
of the stars themselves. This problem is addressed in the 
next subsection where I calculate 95% CL lower limits on 
the expected cluster number density, which in turn provide 
95% CL upper limits on the contribution of VLM stars to 
the halo DM. 

Another important consideration is whether one ex- 
pects HST to be able to resolve individual stars within 



clusters. Since BFGK filter their observations to search 
for only stellar-like sources it is conceivable that a signif- 
icant fraction of the total number of possible VLM star 
sources have been discarded because they reside in the 
cores of unresolvable clusters. Since the number of field 
galaxies far exceeds the number of stars at the typical 
limiting magnitudes of the HST fields, it may prove ex- 
tremely difficult to analyse these fields for the presence of 
unresolved cluster cores. 

Using Eq. 8 as a lower limit on the cluster distance one 
finds a typical angular separation between neighbouring 
stars (as projected along the line of sight) for a cluster 
with radius R of 



-1/2 



£ 0.9/^ 3 /c 1/3 (y 



1/3 / M \ -5/6 



4xl0- 6 sr7 V1O 4 M 



R 



iopcy vo.imq 



1/2 



arcsec, (9) 



compared to a HST pixel resolution of 0.1 arcsec. Since 
one expects clusters to be centrally concentrated one an- 
ticipates Eq. 9 overestimating the angular separation of 
stars near the core and underestimating the separation 
near the cluster edge. Therefore it seems that a significant 
fraction of cluster sources may appear unresolv ed. This 
fraction will be assessed more rigorously in Sect. 3.3 . 

For the cluster scenario to be viable it must be consis- 
tent with existing dynamical constraints on suc h ob jects. 
These constraints, which are reviewed in Sect. 3.4, pro- 



vide upper and lower bounds to both the permitted clus- 
ter mass and radius. The maximum contribution of such 
clusters to the halo is therefore also potentially bounded 
by these constraints. 



Finally, in Sect. 3.5, I show that this scenario is com- 
patible with MACHO microlensing observations, and I 
use the MACHO limit on the contribution of VLM stars 
to compute the allowed ratio of clustered to unclustered 
stars. Potentially, this ratio provides the strongest con- 
straint on the cluster scenario. 

3.2. Space density of VLM clusters 

In Sect. ^| limits were placed on the space density of un- 
clustered VLM stars from HST number counts, assuming 
the stellar distribution obeys Poisson statistics. However, 
since clusters are extended objects, one can not directly 
apply the same assumption to them, since one may have a 
cluster which is only partially within the field of view (in 
fact this will be the case generally for the HST fields due 
to their small solid angle). 

However one can use Poisson statistics to assess the 
likely number of cluster centres within the field of view. 
One can therefore extend this idea by imagining a vol- 
ume centred on the actual field volume, but whose radius 
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of cross section along the length is always larger by an 
amount R (the cluster radius) and whose length is also 
larger by an amount R at either end (i.e. by an amount 
2R overall). Limits on the number of cluster centres ap- 
pearing within this larger volume necessarily correspond 
to limits on the number of clusters which appear either 
wholly or partially within the actual survey field volume. 

The description above applies to a cylindrical survey 
volume (i.e. a circular field of view). The size of the larger 
volume and number-count limits derived from it are sen- 
sitive to the geometry of the field volume on which it is 
centred (i.e. the shape of the field of view), which for the 
HST is 'L'-shaped rather than circular. However, for most 
of the cluster radii considered in this paper their angular 
sizes are always much larger than the angular size of the 
field of view itself, so the sensitivity to the shape of the 
field of view is generally very small, and so for simplicity 
I assume a cylindrical field volume geometry throughout. 

Another important consideration when placing limits 
on cluster number densities is the statistical dependency 
of neighbouring fields. In principle, sufficiently large clus- 
ters could appear in more than one field and so any limits 
would need to take such correlations into account. Even 
if one is placing limits based on the absence of clusters 
within the observed fields one still needs to consider this 
problem. One can see why by considering two neighbour- 
ing fields which are separated by less than a cluster diam- 
eter. The prescription above would oversample the actual 
volume of space required to place limits on the cluster den- 
sity (i.e. there would be a volume of space which would be 
common to the periphery of both fields and would there- 
fore be counted twice), resulting in an overestimate of 
the expected number of clusters. As mentioned in Sect. || 
fields 1 and 2, whilst not overlapping, are nonetheless too 
close for most of the cluster radii considered here and so 
field 1 , which samples a smaller amount of halo mass than 
field 2 (using Eq. ||), is discarded. Similarly, field 19, which 
partially overlaps with field 20, is also discarded leaving 
n = 20 fields. For these remaining fields the spatial separa- 
tion between nearest-neighbours ranges from 150 — 960 pc 
for clusters comprising 0.2-Mq VLM stars (with a me- 
dian separation of 340 pc), and from 50 — 310 pc for clus- 
ters comprising 0.092-Mq stars (with a median value of 
110 pc). The statistical dependency of neighbouring fields 
therefore only becomes important for the very largest clus- 
ter radii considered in this section. I therefore ignore such 
correlations and assume the fields to be statistically inde- 
pendent. 

A lower limit on the expectation number of clusters 
with mass M and radius R, comprising VLM stars of mass 
to* > to, which either wholly or partially appear within 
the n = 20 HST fields is (c.f. Eq. 3) 



N C (M, R,>rn)>^-J2 KsloM m )' 



M 



(10) 



where 



7T I — h R p(x,k,bi)dx. 

'dmin,i(m)-i? V 1 J 



(11) 



0, Mi 



halo.-/ 



Mi 



halo 



in Eq. 



as re- 



Note that for R 
quired. 

If one assumes that the HST fields contain no cluster 
stars at all then the 95% CL upper limit on the allowed 
halo fraction is /h = 3f c N~ x ; 3 being the 95% CL upper 
limit on the true average number of clusters when none are 
seen, assuming statistical independency of neighbouring 
fields and Poisson statistics. Note from Eq. 10 that the up- 



per limit /h is independent of f c since it represents a limit 
on the absolute number of clusters. If f c is significantly 
less than unity stronger constraints on the VLM density 
will come from the limits obtained in Sect. [2] for unclus- 
tercd VLM stars. This issue is considered in more detail in 



Sect. 3.5 where I consider the efficiency with which VLM 



stars need to be clustered. 




Fig. 1. 95% CL upper limits on the allowed halo fraction /h, as 
a function of cluster mass M and radius R, assuming a cluster- 
ing fraction / c = 1. The solid grid shows the limits for clusters 
comprising VLM stars of mass m = 0.2 Mq and the dotted 
grid for stars at the hydrogen-burning limit, m = 0.092 Mq. /h 
is normalised to a local halo DM density po = 0.01 Mq pc~ 3 , 
and fh = 1 is demarcated by the plane on which the skirting 
around the perimeter of the grid is centred. The region of the 
grid at low cluster masses and large radii, which dips below 
this plane, is constrained by HST number counts, whereas the 
region rising above it denotes cluster parameters which can 
provide all of the halo DM whilst remaining compatible with 
HST observations. 



The limits on /j, for a range of cluster masses and radii 
are shown in Fig. [j], assuming f c = 1. The solid grid de- 
picts the limits assuming the clusters comprise VLM stars 
of mass to = 0.2 Mq, where as the dotted line shows 
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the limits for clusters comprising hydrogen-burning limit 
stars (m = 0.092 Mq). The choice of cluster parameters 
is guided by the dynamical considerations discussed in 
Sect. |3.4| . The plane at fh. = 1 is shown to guide the eye 
in seeing which part of the parameter space is constrained 
and which is not. The most strongly constrained region 
occurs for low cluster masses and large radii, though com- 
binations of cluster parameters for which M < 4ttJ>R 3 /3 
(where ~p is the average halo density within the survey 
volumes) are unphysical since they denote clusters whose 
average density is less than that of the local halo back- 
ground. (Such 'clusters' might instead best be regarded as 
voids.) For interesting (i.e. significant) values of f c and fh 
this unphysical region includes the case where the average 
angular separation between neighbouring stars, given by 
Eq. 9, exceeds the HST field of view. Whilst HST provides 
strong limits on clusters of low mass and large radius, it 
is clear from the figure that a significant fraction of the 
parameter space is not constrained by HST observations. 

It is apparent from Fig. [I] that the relationship between 
fh and M is linear, which is simply a consequence of the 
fact that f h ot N- 1 oc M from Eq. [l0|. Whilst the rela- 
tionship between fh and R is not a simple power law one 
can crudely approximate the condition for which fh can 
exceed 1 as 



R<,a 



M 



10 4 M 



© 







(12) 



where for 0.2-Mq stars a ~ 0.2 pc and j3 ~ 1.05, and 
for 0.092-Mq stars a ~ 1.3 pc and j3 ~ 0.75. If instead 
of providing all of the halo DM one merely requires that 
clusters provide at least 40% of the DM, the halo frac- 
tion indicated by the MACHO experiment for a loc al hal o 
normalisation pa = 0.01 Mq pc -3 (Alcock et al. 1997 ), 
then one requires a ~ 0.4 pc and (3 ~ 0.97 for 0.2-Mq 
stars, or a ~ 2.0 pc and (3 ~ 0.73 for hydrogen-burning 
limit stars. These inequalities should be compared to the 
requirement noted above that clusters must represent lo- 
cal density enhancements rather than decrements, which 
requires 



R < 62 



M 



10 4 Me 



© 



1/3 



0.01 Mq pc" 



-1/3 



pc. 



(13) 



The limit should also be compared to the dynamical con- 
straints to be discussed in Sect. 



3.1 



It is important to note that these limits arc fairly in- 
sensitive to the actual number of VLM stars observed 
in the HST fields, provided the observed number is rel- 
atively small (e.g. A^lm < 100). This is because if even 
just a portion of one cluster enters one of the fields of 
view one would generally expect HST to detect hundreds 
or even thousands of stars in that field due to the high 
cluster surface density. For the same reason these limits 
are also relatively insensitive to the surface-density pro- 
file of the cluster, at least within the range of reasonable 



surface-density profiles (i.e. similar to those inferred from 
the surface-brightness profiles of visible clusters). 

3.3. The question of resolvability 

The previous subsection dealt with limits imposed from 
the absence of clusters within the HST fields. However, 
compact clusters can escape detection even if they do ap- 
pear in the fields, provided their surface densities are suf- 
ficiently high as not to allow them to be resolved, or pro- 
vided the total surface density through all clusters along 
the line of sight is sufficiently large. The rigorous point- 
source selection criteria of BFGK means that any unre- 
solved portions of clusters are discarded, in which case 
only some fraction / ros < 1 is detectable. Thus the effective 
number of detectable clusters is N c>e s = (f ICS )N c < N c , 
which implies that the effective allowed halo fraction is 
/h.off oc (f ICS )- 1 N ( T 1 > fh, where (...) denotes averaging 
over fields. Thus if / rcs = for a particular set of clus- 
ter parameters (M, R) then HST point-source counts do 
not place any limits on the allowed halo fraction in such 
clusters. 

The effect on resolvability due to several clusters 
aligned along the line of sight can be estimated from the 
ratio of the average cluster surface density to that between 

^min and G?max- 



7 = /c/h 



7tR 

M 



p(x, i, b) dx. 



(14) 



Assuming the distribution of clusters on the sky is Poisso- 
nian, their sky-covering factor is 1 — exp(— 7). Therefore 
the surface density measured along a line of sight through 
a cluster will be enhanced over that expected for the single 
cluster alone by a factor which is on average 



= 7[1 - exp(-7)]~ 



(15) 



The median line-of-sight surface density between d m i n and 
d max for the 20 HST fields is 25.2 Mq pc" 2 for 0.092-Mq 
stars and 64.6 Mq pc -2 for 0.2-Mq stars, corresponding 
to values for t a of 1.45 and 2.34, respectively, assuming 
10 4 -Mq clusters with a 10-pc radius and f c = fh = 1. 

The resolvable fraction of an individual cluster is sen- 
sitive to its assumed surface-density profile, as well as its 
mass and radius and the distance at which it is expected 
to be observed. I assume here that the cluster surface den- 
sity follows the surface-brightness profile of many observed 
globul ar clusters, which are well described by the King 
( 1962 ) surface-brightness law. Thus the surface number 
density u as a function of cluster-centric radius r is as- 
sumed to be given by 



1 + (r/r c ) 







1 (r < R) 
(r > R) 



(16) 
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where ao = a(r = 0) denotes the central surface density 
and r c is the projected cluster core radius. Since the inte- 
grated mass of the cluster must be M the central surface 
density is 



(M/m) 



a 2 7ri? 2 ln(l + a 



(17) 



where a = r c /R. For the purpose of simplification I adopt 
a — 0.1 for all clusters, a value which is typical for visible 
clusters, though in reality there is a large dispersion about 
this value. 

Clusters appear resolved in the HST fields provided ao 
is less than 



0~Tes (2x$ rcs ) 

4.3 x 10 4 | — I P<- 



5 kpc 



(18) 



where 9 Tea = 0.1 arcsec for the HST. If ao > <J r es then 
one expects that at least some of the cluster is unresolv- 
able. From Eq. |l6| to 18, the mass fraction of a cluster at 
distance x which can be resolved is 



U{x) = 1-- 



min [r ( a — cr rcs / e a ) , R] 



2irra(r) dr 



1 — min < 1 



ln[e a (T /a Ica (x)} 
ln(l + a" 2 ) 



(19) 



if e&cro > a ICS or unity otherwise. Note the factor e CT from 
Eq. |l5|, which enters in the upper integral limit. This is 
because, as argued above, one expects the observed line- 
of-sight density to be e a o, rather than just a. 

Thus, for HST field i, the expectation value for the 
resolvable fraction of clusters is 



(/res, 2) 



X^ max ' 1 x 2 p(x) Ax 



(20) 



where p is given by Eq. 13. Hence, the effective number of 
detectable clusters within the n — 20 fields is given by 

Eq. ^ With MUo,i - (fr e s,i)MUo,i- 

Figure |2| shows the upper limits on the effective halo 
fraction /h, c ff = 3 f c N~ eS , assuming f c = 1. Unlike the 
plot in Fig. [l], the surface does not strictly represent a 
95% CL constraint, since for that one would need to eval- 
uate the joint 95% CL variation in the cluster number 
N c and in the resolvable fraction / rcs . Such an evalua- 
tion would require detailed Monte-Carlo simulations of the 
cluster distribution for each parameter set (m, M, R), and 
for each field i. 

Figure || shows that the problem of cluster resolvabil- 
ity only becomes important for high-mass, compact clus- 
ters, but that it can have a very strong effect on detection 
in this regime. For the highest mass and most compact 
clusters considered the resolvable fraction (/ res ) rapidly 




Fig. 2. Upper limits on the effective allowed halo fraction /h, e ff 
as a result of taking into account the expectation fraction of 
cluster stars which can be resolved (/res). The assumed pa- 
rameters are the same as for Fig. ^ Note that for low cluster 
masses and large radii the allowed fraction is the same as for 
Fig. [I], whilst for large masses and small radii /h, c ff becomes 
asymptotically large as / rcs — > 0. 



approaches zero, so that no point-source observation of a 
VLM star is expected, even if clusters appear within the 
field of view. The density enhancement e a due to the su- 
perposition of clusters is found to have little effect on /h, 
since t a only becomes significantly larger than unity for 
clusters with large radius (in which case e a — > 7 in Eq. ^|) , 
and these clusters anyway have such low surface densities 
that the enhancement has little consequence. 

As for Fig. 0, one can crudely parameterise the portion 
of the plot which permits /h, e ff to be unity by using Eq. ^ 
with a ~ 1.0 pc and (3 ~ 0.73 for 0.2-Mq stars, or a ~ 
1.7 pc and (3 ~ 0.69 for 0.092-Mq stars. If one instead 
simply demands consistency with MACHO (/h, e ff = 0.4) 
then the allowed range is slightly larger with a ~ 1.2 pc 
and /3 ~ 0.74, or a ~ 2.5 pc and (3 ~ 0.70 for O.2-M or 
0.092-Mq stars, respectively. 

3.4- Dynamical considerations 

One can provide strong constraints on the permissible 
mass and radius of clusters which contribute significantly 
to halo DM by considering their dynamical effects on the 



visible stellar population (e.g. Lacey & Ostriker 1985; Carr 
& Lacey |1987| ; Moore |1993| Moore & Silk |1995[ ). Such con- 
siderations lead to the conclusion that viable cluster pa- 
rameters (M, R) are bounded above and belo w, form ing an 
'island' in pa rameter space (Kerins & Carr 1994 ; Moore 
& Silk [l995 ). This is important because it means that 
one cannot invoke arbitrarily massive and compact clus- 
ters in order to force agreement between various observa- 
tions. Furthermore, in principle cluster parameters which 
are compatible with dynamical limits may be incompati- 
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ble with source-count and MACHO limits, thereby ruling 
out the entire scenario. 

The issue of dynamical constraints is somewhat com- 
plicated by the fact that the limits depend (sensitively in 
some instances) on Galactic as well as cluster parameters. 
Whilst I do not explicitly present the dependencies of the 
constraints on these other parameters (these can be found 
in the references cited above), I assume the following val- 
ues for the relevant parameters: a typical Galactocentric 
distance for clusters in the vicinity of the HST fields of 
8 kpc (in fact most of the fields are further from the Galac- 
tic centre, where their dynamical effects would be weaker 
than assumed here); a cluster halo fraction / c /h,eff = 0.4 
(required to provide consistency with MACHO observa- 
tions); an average halo mass density po — 0.01 Mq pc -3 , 
giving a cluster number density f c fh,efipo/M; a halo core 
radius a — 5 kpc; a one-dimensional halo velocity dis- 
persion of 156 km s _1 ; and a Galaxy age of 15 Gyr. The 
dependency on Galactic parameters means that the dy- 
namical limits presented here do not represent firm limits 
in the strictest sense, but should nonetheless help in as- 
sessing the viability of the cluster scenario. (Note that 
some of the authors above express constraints in terms of 
the cluster half-mass radius R1/2 — 0.3 R, assuming the 
surface-density profile of Eq. [l6] with r c = 0.1 R.) 

An upper limit on the radius of clusters comes from 
the requirement that they do not disrupt one another 
due to collisions occurring within the lifetime of the 
Galaxy. This implies R ^ 100 pc for the parameters 
adopted above. Clusters also need to avoid tidal dis- 
ruption due to the differential force from the Galac- 
tic potential acting across the cluster diameter; a con- 
straint which requires R <> 40(M/10 4 Mq) 1 / 3 pc. Since 
at any time there is always a finite fraction of cluster 
members whose velocities exceed the cluster escape ve- 
locity, clusters slowly evaporate over time. That they 
do not evaporate by the present day requires R <; 
0.92 pc (M/10 4 Mq)- 1 / 3 ^^ M )ln(O.4Af/m)] 2 / 3 . 

The final constraint comes from considering the effect 
of close or direct collisions between da rk clu sters an d ob- 
served diffuse globu lar cluster s (Ca rr 197S ; Moore 1993| ; 
Moore & Silk |l995| ). Moore ( |l993[) showed that the re- 
quirement that observed diffuse clusters show no sign of 
disruption today leads to the limit M ^ 10 3 Mq for 
dark clusters with small radius, though this value is de- 
rived from a rather unusually diffuse cluster that may 
have atypical properties which result in too strong a con- 
straint. A more robust limit from other less diffuse and 
more representative clusters demands M & 5 x 10 4 Mq 
for /c/h,eff = 0.4 (B. Moore, private communication). For 
clusters with larger radius Moore & Silk (1995) derive a 



limit corresponding to R > 9 (M/10 4 M Q )^ T pc. 

These limits, for an assumed star mass m — 0.2 Mq, 
are projected onto the plane /h, c ff = 0.4 in Fig. || (bold 
lines), together with the corresponding HST limits from 
Fig. H, which assume a clustering fraction f c = 1. Re- 




Fig. 3. Comparison of HST limits, MACHO observations 
and dynamical constraints, assuming f c — 1 and a VLM 
star mass of 0.2 Mq. The upper limits permissible from HST 
counts (shown in Fig. ^ for which /h, c ff > 0.4 are projected 
onto the plane /h, e fi — 0.4; the central value inferred from 
MACHO LMC observations for a local halo normalisation 
of 0.01 Mq pc -3 . Also projected onto the same plane are 
the dynamical limits discussed in the main text (bold lines), 
which are sensitive to Galactic as well as cluster parameters. 
Note that there exists a small region of parameter space with 
M ~ 4 x 10 4 Mq and R ~ 3 pc which is simultaneously com- 
patible with all of the limits. The overlap is marginally larger 
for lower mass VLM stars. 



gions of parameter space for which HST counts permit 
/h.eff > 0.4 are projected onto the plane /h, c ft = 0.4, cor- 
responding to the halo fraction preferred by MACHO mi- 
crolensing observations towards the LMC (Alcock et al. 



1997). Thus the flat region of the (M, i?) plane at small 



cluster radii and large masses simultaneously satisfies MA- 
CHO observations and HST counts. It is clear that this re- 
gion and the region bounded by the dynamical limits are 
almost mutually exclusive, with only a small area of the 
(M, R) plane satisfying all three requirements. This is due 
to the strong limits inferred from the absence of dynamical 
effects on presently-observed diffuse globular clusters, and 
from evaporation considerations. The surviving parameter 
region is characterised by a cluster mass M ~ 4 x 10 4 Mq 
and radius R ~ 3 pc and is relatively insensitive to the 
VLM star mass within the range considered in this paper. 
Clearly, for the cluster scenario to be viable M and R have 
to be rather finely tuned, though it is interesting that the 
allowed parameter space coincides well with the baryon 
Jeans mass at the recombination epoch, as well as to the 
predictions of some baryonic DM formation theories (e.g. 
Ashman |l990| ). 



3.5. Microlensing signatures 

The effect of clustering halo objects can be to increase 
or decrease the observed microlensing rate, depending on 
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whether or not the observer's line of sight intersects with 
more or less than the expected average number of clusters. 
One can therefore derive firm upper and lower limits on 
the halo lens fraction by assuming it to be clustered. 



Maoz (1994) has shown that there are two key fea- 



tures that should help to discriminate between microlens- 
ing events in the clustered and unclustered cases. The first 
is that the microlensing optical depth through a cluster is 
t ~ MR 2 /mR 2 , where R c = [AGmx(L - x)/c 2 L] 1 / 2 is 
the Einstein radius and L is the distance to the lensed 
source. Since this can be a factor 10 4 higher than in the 
unclustered case for very massive (M ~ 10 6 Mq) com- 
pact [R ~ 1 pc) clusters, Maoz notes that one expects 
to observe angular correlations between events on a scale 
8 ~ R/x. Maoz claims that for clusters of this mass a 
sample as small as 10 events would be sufficient to rule 
out such clusters from providing all the halo DM at the 
95% CL, if such correlations are not observed. Assuming 
that angular correlations are observed, one should also 
expect nearby events to have similar timescales, assuming 
the lenses have similar masses. This is because the events 
would be due to objects residing in the same cluster, so 
they would be at the same distance from the observer and 
would have similar velocities. These conclusions are how- 
ever based on cluster mass scales which are now ruled out 
by the dynamical considerations discussed in the previous 
subsection. The prospects for detecting such signatures 
from clusters w ith M ^ 10 4 Mq are much less optimistic 
(Metcalf & Silk 1996 ), especially if they comprise only 40% 
of the DM. 

Whether clustering could significantly affect the statis- 
tics for the halo DM fraction /h depends upon how many 
clusters one expects to see within the solid angle of ob- 
servation. The MACHO LMC search covers a solid angle 
of 11 deg 2 (Alcock et al. |l996|) , so the expected number 
of clusters between the observer and the LMC (I — 280°, 
b = -33°, L = 50 kpc) is 



N, 



cxp 



where 



j = 



1.3xl0 4 / c / h 



M 



10 4 Mr. 







11 deg' 



Poo," Jo 



x 2 p(x, I, b) dx 



(22) 



for a halo with a density profile given by Eq. [l]. Thus, 
for clusters which are simultaneously compatible with dy- 
namical limits and HST source counts (implying M ~ 
4 x 10 4 M Q from Fig. |) one expects N cxp > 1 300 to be 
within the MACHO field of view if f c = 1 and f h = 0.4. 

If the clusters all have the same mass and comprise 
objects with the same MF, then the observed rate will de- 
pend only on N exp . Assuming the cluster distribution on 
the sky obeys Poisson statistics, the probability of there 
being N clusters within the field of view when one ex- 
pects N cxp (in which case one would infer a halo frac- 



tion fhN/N exp when the true fraction is /h) is simply 
P(N, N cxp ) = (N c N xp /Nl)eM-N cxp ). 

The MACHO data yields a likely halo fraction f h ~ 0.4 
if the lenses are assumed to be unclustered. The in- 
ferred value for fa would be the same for the cluster 
scenario provided N = N cxp , otherwise one should take 
/h ~ 0.4 N cxp /N. The l-er variation on N exp /N is less than 
3% for N cxp = 1 300, so the inferred halo fractions for the 
clustered and unclustered regimes are virtually identical 
in this case. Only for / c Si 0.1 do Poisson fluctuations in 
the number of clusters become important, though in this 
case the microlensing signal is anyway dominated by the 
smoothly distributed VLM stars. The conclusion therefore 
is that for cluster mass scales compatible with dynamical 
limits their microlensing signatures are indistinguishable 
from the unclustered case. 

It is interesting to consider what the inferred halo frac- 
tion from microlensing, together with the HST source- 
count limits, imply for the likely value for the cluster- 
ing efficiency f c . The halo considered in this paper has 
a microlensing optical depth r — 5.6 x 10~ 7 towards 
the LMC, assuming it is completely comprised of lenses. 
The 95% CL lower limit on the me asured optical depth 
is r 95 ~ 1.5 x 10~ 7 (Alcock et al. 1997 ), compared to 
the optical depth contribution expected from all non-halo 
components of r ^ 5 x 10~ 8 . Subtracting the contribution 
from these components gives a lower limit on the halo 
fraction of f h > (1.5 - 0.5)/5.6 = 0.17. From Sect. | the 
maximum contribution from unclustered zero-metallicity 
VLM stars to the halo is 1.4% at the 95% CL. Therefore, 
consistency between MACHO and HST observations re- 
quires that /c > 1 - (0.014/0.17) = 0.92. That is, one 
requires a present-day clustering efficiency of 92% or bet- 
ter. This is certainly a very strong demand for the cluster 
scenario to meet. If one adopts the central value for the 
measured optical depth (r = 2.9 x 10 -7 ) then the required 
clustering efficiency must be at least 97%. 



, (21) 4. Discussion and conclusions 



In this paper I have considered the constraints on the con- 
tribution to the halo dark matter (DM) of a population of 
zero-metallicity, hydrogen-burning stars with mass below 
0.2 M Q (VLM stars). 

Though there are already strong constraints on the 
VLM star scenario, this present work is motivated by 2 
considerations: (1) previous studies have adopted V — I 
colour cuts for their star samples which effectively elimi- 
nate any chance they have of discovering zero-metallicity 
VLM stars; (2) previous studies only consider limits on 
a smoothly distributed population of VLM stars, where 
as some theories of baryonic DM formation predict that 
the stars should be grouped into globular-cluster configu- 
rations. 

Using dat a from 20 of the 22 HST fields obtained by 
Gould et al. ( 1996) ) , together with the photometric pre- 
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dictions of Saumon et al. (1994) for zero-metallicity VLM 
stars, I find that the contribution from a smoothly dis- 
tributed population of such stars to the halo can be no 
more than 1.4% at the 95% confidence level on the basis of 
75 candidate VLM stars with 1.2 < V — I < 1.7 appearing 
in the fields when at least 6 310 are predicted. In reality, 
the true fraction is likely to be less than this value since 
many of the candidates may belong to the spheroid or disc 
components. This limit is comparable to previous analy- 
ses for stars of non-zero metallicity and therefore the in- 
escapable conclusion is that any smoothly distributed pop- 
ulation of VLM stars, regardless of its metallicity, makes 
at best only a tiny contribution to the halo DM. 

Clustering allows the possibility of much larger fluctu- 
ations in the expected number of stars appearing in the 
fields. Additionally, highly compact clusters may not be 
completely resolvable, thereby decreasing the number of 
available point sources. These two effects can permit a 
halo fraction in clusters which is compatible with the halo 
fraction inferred by the MACHO gravitational microlens- 
ing experiment provided that the cluster mass M and ra- 
dius R satisfy the inequality R <> 1.2 (M/10 4 M Q )°- 74 pc. 

One also requires that the cluster scenario satisfies the 
dynamical constraints which exist on the allowed mass and 
radius of clusters. Comparison of the dynamically allowed 
region with the region which satisfies both HST and MA- 
CHO observations reveals that only a very small portion of 
parameter space, characterised by M ~ 4 x 10 4 Mq and 
R ~ 3 pc, satisfies all three requirements. However, the 
implicated cluster mass is intriguingly close to th e mas s 
predicted by some DM formation theories (Ashman |l990 ), 
and to the baryon Jeans mass at the cosmological epoch 
of recombination. 

For a cluster mass M is 10 4 Mq, which is required by 
dynamical arguments, the effect on microlensing statistics 
is indistinguishable from the unclustered case, and in par- 
ticular the inferred halo fraction should be within a few 
percent of the value inferred by assuming the DM dis- 
tribution to be unclustered. However, MACHO and HST 
observations place a strong limit on the efficiency with 
which VLM stars need to be clustered in order to remain 
compatible with both surveys. This limit corresponds to 
a 92% present-day clustering efficiency and therefore pro- 
vides a very stern test of the scenario. 

Whilst recent microlensing results seemingly provide 
strong evidence for a substantial baryonic contribution to 
Galactic halo DM, when taken together with other ob- 
servational and theoretical constraints it has become in- 
creasingly difficult to provide a unique baryonic candidate 
which can simultaneously explain the high microlensing 
fraction and the event timescales. One way out is to in- 
voke a substantial modification in the shape of the halo 
though, for this to work, the microlensing results require 
a component resembling something closer to a maximal 
disc, leaving the status of any halo (and the role of non- 
baryonic DM on Galactic scales) much reduced. Another 



option is to attribute the lensing events to some non- 
baryonic candidate, such as primordial black holes, though 
in this case one requires an additional non-baryonic can- 
didate to explain the rest of the halo DM. Here I have 
shown that it is still possible, if only barely, to construct 
a baryonic scenario which is compatible with all known 
constraints and which does not require a major modifi- 
cation in the halo dynamics. It is also conceivable that 
the scenario could provide an explanation for the recent 
detections of faint extended emission around the edge-on 
spi ral ga laxy NGC 5907 (Sackett et al. 1994 ; Lequeux at 
al. 1996|) which, under the assumption of constant mass- 
to-light ratio, appears to trace the distribution of a halo. 
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